Far field imaging by a planar lens: diffraction versus superresolution 
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We resolve the long standing controversy regarding the imaging by a planar lens made of left- 
handed media and demonstrate theoretically that its far field image has a fundamentally different 
origin depending on the relationship between losses inside the lens and the wavelength of the light 
A. At small enough A the image is always governed by diffraction theory, and the resolution is 
independent of the absorption if both Ime ^ 1 and Im^ <C 1. For any finite A, however, a critical 
absorption exists below which the superresolution regime takes place, though this absorption is 
extremely low and can hardly be achieved. We demonstrate that the transition between diffraction 
limited and superresolution regimes is governed by the universal parameter combining absorption, 
wavelength, and lens thickness. Finally, we show that this parameter is related to the resonant 
excitation of the surface plasma waves. 

PACS numbers: 42.30.-d,73.20.Mf,42.25.Fx 
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The left-handed medium (LHM), introduced by 
Veselagoi in 1967 as a medium with simultaneously neg- 
ative and real ji and e provides a negative refraction at its 
interface with a regular medium (RM). This effect allows 
creation of a unique imaging device, sometimes called 
the "Veselago lens", formed by a planar slab of LHM, 
with its refractive index and impedance ideally matched 
to the surrounding RM (Fig.l). Interest in the planar 
lens significantly increased after the work by PendryS, 
who suggested that the planar lens can in principle fo- 
cus all Fourier components of a 2D image, introducing 
a term "the perfect lens" . However, taken literally, this 
statement contradicts Electrodynamics: since there is no 
source at the focal point the field of the source cannot 
be exactly reproduced in the vicinity of the focal point. 
As it was pointed out in Rcfs.[3A_5,], in the absorption- 
less limit (Ime — Im/i = 0), the solution proposed ir^ 
exponentially diverges inside a 3D domain between two 
foci (see Fig.l), and therefore it cannot be a solution of 
Maxwell's equations. Pokrovsky and Efros^ proposed a 
diffraction theory that should be valid at large fco = w/c 
and that does not contain any superresolution. Haldane^ 
explained that the superresolution should be connected 
to the resonance of the plasma waves on both sides of 
the slab, but he added that the theory should be regular- 
ized. The regularization due to a very small absorption 
inside the LHM, and its implication for resolution were 
analyzed by many aut hor s^i^i^ i i o d i d ^ _ Near field behav- 
ior has been studied since 1994i^, both theoretically 
and experimentallyi^ii^iiiii^. The goal of this work is to 
develop a quantitative criterion of applicability of diffrac- 
tion theories for planar lens systems. 

It has been shown that the far-field imaging of the 
planar lens has the following properties: 

(a) The stationary solution of the Maxwell equations 
does not exist when lnie=lnifi—0. When absorption is 
small, but non-zero, there are two "resonant" regions 
shown in Fig. 1; the fields inside resonant regions di- 
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FIG. 1: (Color online)Planar LHM-based lens with two real 
foci outside and inside the slab. S is the source, a is the 
distance from the source to the slab, d is the width of the 
slab. Dark (Red) regions with dashed boundaries are the 
resonant regions. The fields diverge in these regions when the 
absorption tends to zero. 



verge as absorption tends to zero. Both foci arc at the 
boundaries of the resonant regions. Therefore the focus 
outside the lens is quasi virtual: an observer to the right 
of the focal point can see subwavelength focus while an 
observer to the left of the focus can see only large and 
highly oscillating fields - an indication of resonantly ex- 
cited surface plasma waves propagating at the back in- 
terface of the lens. The existence of such a quasi-focus 
does not contradict general theorems. 

(b) Milton et al^ showed that total absorption of elec- 
tric energy inside the slab oc (Im e)(2a/d-i)^ where a is 
the distance from the source to the lens and d is the width 
of the lens. If a < d/2, the resonant regions overlap. In 
this case the total absorption inside the lens tends to in- 
finity as Ime tends to zero. Moreover, the electrostatic 
lens cloaks two dimensional dipoles rather than imaging 
them "perfectly" (See also Ref.[ll]). 

On the other hand, 2D and 3D diffraction theories 
demonstrate a real focus with different widths in the 
lateral and in the longitudinal directions. According to 
these theories, fields near the focus are the universal func- 
tion of r/A, where r is the radius- vector with an ori- 
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gin at the focal point^i^'^. Although diffraction theory 
uses spherical waves rather than plane waves, results of 
both 2D and 3D diffraction theories coincide with so- 
lutions in the plane wave representation with evanescent 
waves omitted. The diffraction theory is known to be the 
first approximation at small deviations from geometrical 
opticE^^-^^. However in the case of the planar lens we 
have a paradoxical situation. If Ime <^ 1 and Imfi <^ 1, 
the width of the focus, as obtained by the diffraction 
theory, is independent of absorption, while the exact so- 
lution has a singularity at zero absorption even at small 
wavelengths. Therefore, the criterion of applicability of 
the diffraction theory cannot be related to a trivial ratio 
of wavelength to width of the slab, and must also depend 
on absorption. 

Podolskiy and Narimanov^ studied the width of the 
focus of the 2D planar lens. They have found that it is 
close to the diffraction limit almost at all cases where 
the wavelength is smaller than the width d of the LHM 
slab. They also made an important observation that the 
boundaries of the diffraction regions depend on the ab- 
sorption. 

First we describe our results. We consider TM po- 
larization with a 2D Green's function as a source. The 
magnetic field with an amplitude h near the source at 
the origin has a form Gr = {i / 4)hH'^\pkQ) , with H 
being the Hankel function h'^^ = Jq -\- iNq, and p — 



y^. This function can be represented in the form 



Hev, where 



ih r'^o expz 
Hp — — 

TT 



(ky + z^/k^ 



k'^ -LUt 



dk (1) 



contains only propagating waves while 
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exp (iky — z^Jk"^ — fcg — iujtj 



|fe|>feo 



dk (2) 



contains only evanescent waves (EW). 

The 2-dimensional diffraction theory22i2^ shows that 
the field near the focus contains only propagating modes. 
If the slab is a,t a < z < d-\- a and Ime <C 1 and Im/it <C 1 , 
the field has a form 



Hliy,z') 



ih f''° exp i [ky + z'y^kg — k"^ — Lut^ 
TT J-ko ^Jkl - k^ 



dk. 



where z' = z — {2d — a) is a distance from the focus in z- 
direction. At z' = one gets H^{y,0) = ihJo{2Try/X) 
Note that the halfwidth of the first maximum of the 
square of the Bessel function is Ay w 0.3A. This small 
value of the width in the lateral direction has been mis- 
interpreted as superlensingSii^. 

As a criterion of transition from the supersolution 
regime to the diffraction regime we choose the ratio 
P = Hev/Hf{0, 0) = Hev/ih, where both magnetic fields 
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FIG. 2; (color online) Dependence of the parameter \P\ 
on normalized lens thickness for S — 10~^ (solid triangles), 
S = 10"^ (boxes), S = 10"" (stars), and S = 10~^* (empty 
triangles); symbols correspond to results of numerical inte- 
grations; dashed and solid lines correspond to approximate 
Eq.® and Eq.® respectively; note that quasistatic Eq.([5)) 
adequately describes the far-field behavior of the system in 
the limit of vanishingly small absorption 



of EW's and magnetic fields given by the diffraction the- 
ory are taken at the focal point z = 2d — a. Then if P 
is small, we have the diffraction regime, while at large P 
we have the regime of supersolution. 

One should keep in mind that the creation of the 
LHM is a difficult and controversial problem. Usually 
the LHMs are metamaterials with periodic or quasi- 
periodic structure^^ and their magnetic permeability /i 
may exhibit strong spatial dispersion [/i = p,{uj, k)^^i^. 
All these problems are outside the scope of this paper. 
Here we consider the basic model - a homogeneous and 
isotropic "hypothetical" slab that has 



e — p, — —I + i6, 



(4) 



where S 1; e = ji — 1 outside the slab. Our goal is to 
show what one can expect in the best case scenario. 
Our main results are: 

(i) The parameter responsible for the far-field (fcod S> 
1) transition between diffraction limited and superreso- 
lution regimes is not described by the wavelength alone, 
but also depends on absorption: 



S = fcorfVlm(eT/x)72 = kody/S. 



(5) 



(ii) At 5 » 1 the system is always described by diffrac- 
tion theory, and 

p^g^ ^ 1 4(iisin5-cos5')exp(-5') ^ ^ 



kod 



(iii) At S* < 1 (and kod > 1) the system may exhibit 
superresolution, with 
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P(5) = -sinh-i(y//(2M)), 



(7) 
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where y/ is the root of equation 
45"* exp y / 



1. 



(8) 



(iv) Finally, when - ln(^/4) > k^d [S <C 1) the system 
exhibits superresolution with quasistatic-like behavior 



P{S) 



■In 



2 ln(^/4) 



cxlnln(l/S'). 



(9) 



P[S) ^ 00 as 5* ^ 0. However, divergence of P{S) is ex- 
tremely slow, and in practice |-P(>S')| ^ 1 is unachievable 
in realistic far-field structures. 
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Note that the limit A — > (fco — > oo) at a fixed 5 and 
the limit (5 — > at fixed fco(A) do not commute. Indeed, 
at small absorption the diffraction regime exists at large 
enough k and that at large k the superresolution regime 
exists at small enough absorption; as described above, 
the solution at (5 = 5 = does not exist. Below we pro- 
vide the derivation of Eas. (|5l6l7l9p . and demonstrate the 
connection between the superresolution and interaction 
of plasnionic surface waves, which disappears at S* > 1. 



To derive Eqs. ()l|5l9p we use transmission coefficients 
for EW's as calculated by Podolskiy and Narimanov^ and 
take into account that iJ^(0,0) = ih. Then 



P = 

D = 



ex.p{—dy/ k'^ — kl)dk 



™ Jko y/P - /cq sinh {y/k'^ - k^e^d) + cosh (^P - k^efid)] 
kie{e + fi) ^ k^l + e^) 



2ey/k^ - k^y/k^ ~ k^Cfi 

At k^d 3> 1 the values of k, that are important in this integral are very close to feg. Namely, 

k^ -kl^ l/d^. 

Then (fc - fco)/fco ^ l/d'^kl < 1, fc^ - fcg _ 2ko{k - /cq), and = e^i I - 2iS [see Eq.g])]. Thus, 

2 r exp{-y)dy 



TTkodJo ^H-iS'V2y2 sinh ^/y^ -f iS'2/2 + cosh y/y^+iS'^/2' 

I 



(10) 



(11) 
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As it can be explicitly verified, in the regime 5^1, the 
denominator of this integral can be further simplified: 
y/TTiSy2^smh + + cosh y/y^ + iS^/2 ~ 

v/iS^TV sinh y/iS^/2, leading to Eq.©. 

If 5 = the integral in Eq. p^ diverges which means 
that there is no solution without absorption. To consider 
the case of small S it is more convenient to start with the 
transmission given by Eq.(2) of Refi^ that is written for 
the case of small abso rption. Introducing a new variable 
of integration y — 2y/k'^ — k^d one gets 



P ^ 



where 
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dy 



TTfcod Jo V(y/(2M)2 + l(l + </)2 exp(j/)) 



, (13) 



(14) 



In the limit of extremely small absorption, this integral is 
dominated by the first term in Eq. , yielding a quasi- 
electrostatic- like regime [Eq.®], predicted by Nicorovici 
et alM; where the results are independent of n Note that 
this quasi- electrostatic regime may exist even at k^d^ 1 
if absorption is small enough. 



For somewhat larger absorptions, the far field regime 
is described by the second term in Eg. lfH)) . leading to 



P = 



1 



dy 



Trkod Jo ^{y/{2kod)^ + l{l + 4S^ey/y^) 



(15) 



At 5 ^ 1 one gets an adequate approximation for the 
integral Eg p^ assuming that 1/[1 + 45^ exp(y)/y^] — 
'd{yf — y), where 19(2;) is the Heaviside Step function and 
yj is given by Eq.®. This leads to Eq.©. 

Fig [2] shows the comparison of our analytical results for 
P{S) for S* <C 1 and S* 1 to the results of numerical 
integration of Eq. (fTO| . 

We now explain the physical meaning of parameter 
5. The condition kod » 1 is not sufficient for the ap- 
plicability of diffraction theory in the planar lens be- 
cause of the resonance interaction of the surface plasmon 
waves'^. These waves exist in TM polarization under the 
condition^i2iiS^ that the ratio k^/e changes sign at the 
interface of the vacuum and the LHM. At (5 ^ 1 it reads 



k'^-kl^ Jk"^ -kl + kl2i5. 



(16) 



The mismatch responsible for the breakdown of resonant 
excitation of surface modes can be related to the term 
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6n. 




FIG. 3: (color online) The normalized lens thickness as a 
function of absorption for constant values of parameter |P|; 
filled triangles, squares, starts, diamonds, and empty triangles 
correspond to |P| = 1, \P\ = 1/2, \P\ = 1/4, |P| = 1/8, and 
|P| = 1/16 respectively. In the part of the plane above the 
curve with a given \P\, the relative contribution of evanescent 
waves is less than |P|. 

kl2iS. At S* ^ 1 one can use Eq. pdj) to find that the 
absolute value of the ratio of the mismatch term to — 
fcp is S*^ 1. Thus, at 5* 1 there are no traces of 
the resonance and regular diffraction theory should be 
appUcable. At S' < 1 it follows from Eq. that 
k'^ — kg ~ y'j /(P. fgnoring in this estimate the logarithmic 
factor in and assuming ~ 1 we resolve that the 
relative mismatch in Eq. (|16p is also of the order of S^. 
Therefore, the violation due to the absorption is small at 
S -^1 and the main features of the resonance should be 
preserved. 



Finally, we use the developed formalism to analyze the 
range of parameters where one could expect superresolu- 
tion (and quasi virtual focus with subwavelength thick- 
ness). These results are summarized in Figl3] which shows 
the dependence of normalized lens thickness d/X as a 
function of absorption S for a set of fixed values of the 
parameter |P|, calculated using direct numerical calcu- 
lation of the integral in Eq. pH)) . One can see that the 
contribution of the evanescent waves at the focal point is 
practically negligible at d/A > 2 at any reasonable value 
of absorption. 

In conclusion, we have developed an approach to cal- 
culate a quantitative measure of superlensing in a planar 
hypothetical LHM-based lens, and used the developed 
formalism to separate the regions of the superresolution 
and diffraction in the far field regime. We demonstrated 
that the limits of absorption S ^ and wavelength 
A ^ do not commute; the former limit yields superlens- 
ing, while the latter leads to diffraction limited behavior, 
which typically dominates the far-field image of realistic 
planar lenses. We demonstrated that if A and 5 are both 
finite, the behavior of the planar lens is described by a 
universal parameter S which depends on both geometrical 
sizes and absorption^ and found analytically asymptotical 
behaviors for S* — s- and 5 ^ oo. A connection between 
the value of S and the existence of resonant excitation 
of plasmonic waves has also been demonstrated. Under- 
standing the onset of the diffraction limit presented in 
our work is important for the further development and 
design of imaging systems with negative refraction. 
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